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INVESTIGATION  INTO  THE  EFFECTS  OF  WEAPON  SETBACK 
ON  VARIOUS  MATERIALS  AND  GEOMETRIES 


I.  INTRODUCTION. 

The  design  of  chemical  munitions  involves  several  factors  and  among  them  is  the 
structural  behavior  of  the  components  when  subjected  to  various  loads.  In  some  cases  the 
requirement  is  that  the  component  withstand  the  maximum  stress  or  weapon  setback  force  while 
others  require  failure  at  minimum  setbacks  yet  survive  the  logistical  cycle.  In  either  event,  the 
structural  integrity  of  the  component  must  be  evaluated  as  part  of  the  development  of  an  item. 

At  present,  this  evaluation  involves  the  calculation  of  the  component’s  structural 
integrity  using  statically  obtained  strength  values,  followed  by  bench  testing,  if  possible,  and  the 
testing  of  munitions.  The  emphasis  is  on  system  testing  which  yields  specific  solutions  for 
specific  configurations,  and  if  the  configuration  is  changed  the  system  must  be  retested.  This 
method  has  resulted  in  increased  project  costs  and,  at  times,  in  schedule  slippages. 

What  is  required  is  a reliable  method  for  analyzing  the  effects  of  the  launch 
environment  on  munition  components.  Therefore,  an  investigation  was  undertaken  to  develop 
such  a method  and,  thereby,  reduce  the  amount  of  munition  flight  testing  and/or  dynamic  bench 
testing  required  with  its  associated  costs. 

This  work  was  divided  into  two  phases,  the  first  of  which  was  to  develop  a 
mathematical  model  which  predicts  the  mechanical  behavior  of  materials  subjected  to  setback. 
This  involved  a literature  search  followed  by  a mathematical  modeling  effort.  The  second  phase 
involved  comparing  the  model  prediction  with  experimental  data  which  were  generated  as  part  of 
this  effort. 

II.  BACKGROUND. 

Plastic  deformation  of  solids  was  first  considered  in  1904  by  Hopkinson  who 
observed  that  iron  and  copper  wires  subjected  to  rapidly  applied  tensile  stress  could  be  stressed 
beyond  the  static  elastic  limit  and  breaking  loads  and  still  remain  in  the  elastic  range,  providing 
that  the  duration  of  such  stress  was  of  the  order  of  0.001  second  or  less.*  But  it  was  not  until 
1941  when  Von  Karman  and  Taylor  each  independently  established  the  theory  of  plastic  strain 
propagation  in  metals  that  this  phenomenon  was  seriously  studied.-  This  theory  assumes  that 
stress  is  a unique  function  of  strain  and  that  this  functional  relationship  can  be  obtained  by 
solving  the  equation  of  motion  for  whatever  case  is  under  consideration.  The  original  paper  used 
the  simple  case  of  a slender  rod  subjected  to  longitudinal  impact.  Two  important  features  of  this 
theory  are  that  the  velocity  of  propagation  is  less  than  that  of  an  elastic  strain  and  that  an  upper 
or  critical  impact  velocity  exists  above  which  the  material  fails. 

Since  the  publication  of  this  theory,  a body  of  research  data  has  grown;  in  general  it 
can  be  placed  into  four  categories: 
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(1)  The  determination  of  the  dynamic  properties  of  various  materials.  The 
materials  tested  have  been  mostly  high-strength  ferrous  alloys3  and  lightweight  alloys  for  use  in 
the  aerospace  industry4  with  a small  amount  of  data  on  annealed  nonferrous  alloys.  In  general,  a 
material's  strength  increases  when  subjected  to  impact  loads  (relative  to  statically  applied  loads), 
but  by  varying  relative  amounts  for  various  materials.  For  example,  cold-rolled  low-carbon  steel 
has  an  ultimate  strength  of  87,750  psi  under  static  conditions  and  1 25,250  psi  under  dynamic 
conditions  or  a 43%  increase  while  24-ST  aluminum  has  an  ultimate  strength  of  65,150  psi  under 
static  conditions  and  68,600  psi  under  dynamic  conditions  or  only  a 5%  increase.  It  was  also 
noted  that  the  strain  rate  affects  both  the  strength  and  ultimate  elongation  of  a material. 

(2)  Various  investigations  into  the  behavior  of  materials  subjected  to  impact  loads. 
This  work  usually  involves  solving  the  equation  of  motion  for  a geometry  of  interest,  having  one 
or  more  “dynamic”  parameters,  usually  viscous  or  flow  parameters.  This  is  followed  by  dynamic 
bench  tests  to  determine  the  specific  value  of  parameters  for  a specific  material. 

(3)  High-energy-forming  technology  for  manufacturing  applications.  This  work  is 
concentrated  in  two  areas;  the  development  of  techniques  for  applying  high-velocity  stress  waves 
and  the  development  of  process  parameters.  The  techniques  developed  include  both  high-  and 
low-explosive-forming,  exploding-wire  (capacitor  rapidly  discharging  through  a fine  wire  which 
vaporizes)  forming,  magnetic-pulse  (quickly  changing  electromagnetic  field)  forming  and 
mechanical-pneumatic  (air-activated-ram)  forming  techniques.5  The  parametric  process  studies 
involved  relating  the  depth  of  draw,  material  thickness,  die  diameter  and  standoff  to  the 
material’s  mechanical  properties  and  the  amount  of  forming  energy  involved.4 

(4)  Physical  descriptions  of  the  effect  of  impact  loading  of  metals.  The 
displacement  mechanisms  for  metal  subjected  to  impact  loads  are  similar  to  that  for  statically 
applied  loads,  that  is,  grain  distortion,  slip*  and  twinning.**  But  in  the  case  of  impact  loads,  the 
deformation  is  greatest  near  the  point  of  application  and  reduces  with  distance  from  that  point.6 
There  is  also  a time  delay  between  the  application  of  an  impact  load  and  the  subsequent 
yielding.  This  phenomenon  has  been  explained  in  part  by  the  presence  of  foreign  atoms  in  a 
metal  lattice  and  grain  boundaries,  both  of  which  constitute  obstacles  to  the  flow  of  metals.1 
Another  obstacle  to  flow  is  strain  hardening  which  is  due  to  the  buildup  of  dislocation  dipolesf 
into  a network  which  anchors  dislocations.7 

In  summary:  (l)The  displacement  mechanisms  for  metals  subjected  to  impact  loads 
are  similar  to  those  for  statically  applied  loads,  but  with  different  distributions,  (2)  there  exists  a 
time  delay  between  the  application  of  impact  loads  and  subsequent  plastic  yielding,  (3)  materials’ 
strength  increases  when  subjected  to  impact  loads,  but  by  varying  relative  amounts  for  various 
materials,  and  (4)  previous  investigations  showed  that  flow  or  viscous  parameters  are  important. 


* Slip  is  the  displacement  of  one  part  of  a crystal  relative  to  another. 

**  Twinning  is  the  sliding  of  one  plane  of  atoms  over  the  adjacent  plane. 

t Dislocation  dipole  is  the  junction  of  one  dislocation  with  another  dislocation  on  an  intersecting  plane. 
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III.  MODELING  EFFECT. 


Based  on  the  information  obtained  during  the  literature  search  a mathematical 
model  lias  been  developed  to  predict  the  time-dependent  behavior  of  materials.  This  model 
assumes  that  materials  behave  in  an  elastic  manner  below  the  proportional  limit  (arbitrarily  set  at 
2%  elongation)  and  in  a viscoelastic  manner  above  the  proportional  limit.  In  addition,  the 
viscoelastic  element  must  have  a time  delay  or,  in  this  case,  a viscous  resistance  feature  which  is 
known  as  a Voigt  element.  The  mechanical  analog  of  the  model  is  given  in  figure  1,  the 
equations  for  which  are  given  to  the  right  of  the  figure. 


a 


a 


Ee 


de 

v7t  + Ee 


Figure  1 . Mechanical  Analog  of  Elastic-Viscoelastic  Model 

The  response  equation  for  this  model  is  o = (2E)e  + tj  * j^g  time-dependent 
strain  for  this  system  is  found  by  solving  this  equation  for  the  model  when  excited  by  the  forces 
due  to  the  launch  environments  (which  can  be  obtained  from  acceleration/pressure  histories  of 
the  weapon  system).  Once  the  strain  relationship  is  known,  the  other  dynamic  values  can  be 
determined.  This  involves  taking  the  Laplace  Transform  of  the  dynamic  equation,  rearrangement 
and  taking  the  inverse  transform  to  find  the  time-dependent  strain.  The  “dynamic  modulus"  is 
then  the  stress  divided  by  the  strain  which  is  constantly  changing  with  time. 


* See  the  glossary  for  definition  of  variables  used  here  and  in  succeeding  sections  of  this  report. 
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The  elastic  and  viscous  constants  of  a material  may  be  obtained  from  a statically 
obtained  true  stress-strain  diagram  by  assuming  that  the  material  behaves  elastically  below  the  pro- 
portional limit  and  in  a viscous  manner  above  it.  The  graphical  representation  of  this  is  shown  in 
figure  2. 


Figure  2.  Idealized  True  Stress-Strain  Diagram 
IV.  EXPERIMENTAL  METHODS. 

To  compare  the  model  predictions  with  experimental  data  a series  of  dynamic  bench 
tests  was  performed  involving  two  geometries  (the  thin  flat  plate  as  shown  in  figure  3 and  the 
solid  pin  as  shown  in  figure  4)  and  four  materials  (AISI  1018  and  1020  low  carbon  steel, 
type  316  austenitic  stainless  steel  and  6061-T6  aluminum  alloy).  In  addition,  the  stress-strain 
diagrams  for  all  four  materials  were  generated  to  provide  elastic  and  viscous  material  constants. 

The  bench  test  machine  used  in  this  work  was  a Veripulse  VP  400  Shock  Machine 
which  produces  a controlled  and  reproducible  acceleration  at  a constant  pulse  rate.  The 
procedure  for  this  testing  involves  securing  a component  holding  fixture  and  weight  to  the  bed  of 
the  shock  machine  and  varying  the  bed  drop  height  to  achieve  the  proper  acceleration  while  the 
neoprene  compression  pads  control  the  pulse  rate.  The  calibration  chart  for  both  pads  is  given  in 
figure  5,  and  the  estimated  standard  deviation  for  each  pad  is  given  in  the  results  section  of  this 
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PLATE 

Figure  3.  Flat  Plate  Configuration 


Figure  5.  Calibration  Chart  for  Shock  Mixture 
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report.  The  acceleration  history  of  this  machine  is  a half  sine  wave  pulse  (aft)  = ap  sin(-^-)t).  The 

“dynamic  modulus”  for  the  various  materials  can  be  estimated  by  solving  the  dynamic  equation 

for  the  l lastic-Viscoelastic  model  when  excited  by  the  wave  pulse,  the  solution  of  which  is  given 

in  appendix  A.  Note  that  for  loads  of  short  duration,  the  modulus  is  dependent  on  the  viscous 

constant  rather  than  both  material  constants  (elastic  and  viscous).  The  stress  on  a flat  plate  can 

then  be  found  by  solving  the  equation  of  motion  for  the  flat  plate  given  in  appendix  B.  The 

stress  on  a solid  pin  can  also  be  found  by  solving  the  equation  of  motion  for  that  geometry,  the 

solution  of  which  is  analogous  to  that  for  a flat  plate,  but  with  different  boundary 

conditions  (a  = — and  p = ~ ^ 

V 2A  2A  )' 

V.  RESULTS.  , 

The  results  of  the  dynamic  bench  test  effort  are  given  in  table  1.  This  work  involved 
five  material-geometry  combinations  and  because  of  the  go-no-go  nature  of  the  test  setup  only 
upper  and  lower  strength  values  are  reported.  The  Elastic-Viscoelastic-model  predictions  for  all 
five  combinations  are  also  given  in  table  1.  A sample  calculation  for  1018-flat  plate  combination 
is  given  in  appendix  C.  For  comparative  purposes,  strength  predictions  based  on  statically 
obtained  material  properties  using  standard  formulas  are  included  in  table  1 . 


Table  1.  Comparison  Between  Measured  and  Predicted  Strengths 


Material 

Geometry 

Bench-tested 

value 

Model 

prediction 

“Static 

prediction” 

1018 

Flat  plate 

647-693  psi 

697  psi 

27  psi 

Carbon  steel 

(.0153  inch  thick) 

1020 

Solid  pin 

1,048-1,082  lb 

1,111  lb 

388  lb 

Carbon  steel 

(.125  inch  diameter) 

6061-T6 

Flat  plate 

508—555  psi 

486  psi 

59  psi 

Aluminum 

(.025  inch  thick) 

6061-T6 

Solid  pin 

1,048-1,1 19  lb 

1,021  lb 

7141b 

Aluminum 

(.1875  inch  diameter) 

316 

Solid  pin 

2,484-2,552  lb 

9,232  lb 

1,8981b 

Austenitic  steel 

(.1875  inch  diameter) 

The  statically  obtained  true  stress-strain  diagrams  for  all  four  materials  are  given  in 
figure  6.  From  these  diagrams  the  dynamic  modulus  (0  = 7-  17)  and  velocity  of  wave  propagation 

" / /n  \ . i ^ I 


[-©1 


were  calculated  and  are  given  in  table  2. 
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Figure  6.  Statically  Obtained  True  Stress-Strain  Diagram 
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Tabic  2.  Calculated  Dynamic  Material  Constants 


Material 

Geometry 

Dynamic  modulus 

Velocity  of  wave 
propagation 

1018  Carbon  steel 

Rat  plate 

psi 

4,358,000,000 

n 

21,367 

1020  Carbon  steel 

Solid  pin 

4,222,000,000 

21,913 

6061-T6  Aluminum 

Flat  plate 

1,632,000,000 

22,840 

6061-T6  Aluminum 

Solid  pin 

U05.000.000 

20,430 

316  Austenitic  steel 

Solid  pin 

27,390,000,000 

54,400 

To  estimate  the  precision  of  the  Veripulse  VP  400  Shock  Machine,  a series  of 
instrumented  tests  was  performed  on  both  sets  of  neoprene  compression  pads.  For 
pad  31-200-110-200,  used  in  the  1020  steel  and  6061-T6  aluminum  solid  pin  tests,  the  standard 
deviation  was  2.9  g acceleration  or  49.01b  load  (F  = ma  = 16.9  lb  X 2.9g  = 49  lb).  For  pad 
31-200-085-350,  used  in  the  flat  plate  and  the  316  austenitic  steel  solid  pin  tests,  the  standard 
deviation  was  4.6  g acceleration  or  21.2  psi  for  the  plates  and  77.7  lb  load  for  the  pin. 

VI.  DISCUSSION  OF  RESULTS. 

For  the  low  carbon  steel  the  Elastic-Viscoelastic  model  predicted  values  slightly 
higher  than  the  upper  bench  test  value  for  both  geometries  tested  (697  psi  versus  693  psi  for  flat 
plate  and  1,111  lb  versus  1,082  lbs  for  the  solid  pin).  For  the  6061-T6  aluminum  alloy  the 
Elastic-Viscoelastic  model  predicted  values  slightly  lower  than  the  lower  bench  test  value  for 
both  geometries  tested  (486  psi  versus  508  psi  for  the  flat  plate  and  1,021  lb  versus  1,048  lb  for 
the  solid  pin).  In  each  case  the  difference  between  predicted  and  measured  values  are 
approximately  the  estimated  standard  deviation  of  the  shock  machine  and  consequently  the 
differences  may  not  be  significant.  In  comparing  the  two  methods  of  predicting  strength  values, 
it  is  seen  that  the  predictions  based  on  the  Elastic-Viscoelastic  model  are  far  superior  to  those 
made  using  standard  material  strength  values  and  formulas. 

For  the  austenitic  stainless  steel  type  316  the  Elastic-Viscoelastic  model  predicted 
values  considerably  higher  than  the  upper  bench  test  value  for  the  solid  pin  geometry  (9,232  lb 
versus  2,522  lb).  The  most  obvious  explanation  for  this  discrepancy  is  that  the  model  as 
proposed  does  not  adequately  account  for  the  behavior  of  316  while  it  does  so  for  low-carbon 
steel  and  the  aluminum  alloy. 
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In  looking  into  the  structure  of  the  three  metals  some  differences  were  found.  Both 
the  low-carbon  steel  and  6061-T6  aluminum  alloy  are  basically  pure  metals  doped  with  small 
percentages  of  other  atoms  to  achieve  certain  desirable  mechanical  properties.  Furthermore,  both 
materials  are  subject  to  strain  hardening.®  Both  the  presence  of  impurities  in  the  crystal  lattice 
and  the  susceptibility  to  strain  hardening  tend  to  produce  barriers  to  metal  flow  which  can  be 
modeled  as  a viscous  drag  element. 


On  the  other  hand,  austenitic  stainless  steel  type  316  is  a mixture  of  iron(Fe), 
chromium  (Cr)  and  nickel  (Ni)  where  the  Cr  and  Ni  atoms  substitute  for  Fe  in  the  crystal  lattice. 
Because  of  the  high  nickel  content  (10  to  14%)  this  steel  does  not  appreciably  strain-harden  and 
is  sometimes  referred  to  as  free  spinning  steel.^  This  relative  lack  of  strain  hardening  has,  in  part, 
been  explained  by  an  “easy  glide”  mechanism  where  the  density  of  dislocations  rises  linearly 
with  plastic  strain.7  In  fact,  the  presence  of  grain  boundaries  and  impurities  would  be  the  only 
serious  barrier  to  metal  flow.  Since  the  austenitic  stainless  steel  does  not  follow  the 
Elastic-Viscoelastic  model  previously  stated,  a new  model  for  this  material  must  be  postulated. 
Based  on  the  above  discussion  concerning  316,  a three-element  model  shown  in  figure  7 is 
proposed.  This  model  involves  an  elastic  element  to  account  for  the  behavior  below  the 
proportional  limit  followed  by  a Voigt  element  to  account  for  grain  boundary  behavior  followed 
by  a friction  element  to  account  for  the  “easy  glide”  phenomenon.  The  equation  of  motion  for 
de 

this  system  is  a = (2E  + 3)e  + rj  whose  equation  is  analogous  to  that  for  the  Elastic-Viscoelastic 
model,  i.e.,  4>~  ~ V' 


Figure  7.  Mechanical  Analog  for  Elastic-Viscoelastic-Plastic  Model 


16 


In  analyzing  the  true  stress-strain  diagram  it  is  seen  that  316  is  clastic  from 
0 to  .01  strain,  viscoelastic  from  .01  to  .039  and  plastic  from  .039  to  59.  For  this  case  the  viscous 

parameter  t?  = ---  00333'^"  ~ ( 1 + .039)  = 3.4  X 10^;  likewise,  <f>  = 5.91  X 10**  and 

c = 13,040  ft/sec.  Substituting  this  value  of  c into  the  equation  of  motion  for  the  solid  pin  yields  a 
predicted  value  of  2,861  lb.  In  comparing  the  three  methods  of  predicting  strength  values  it  is 
seen  that  the  prediction  based  on  the  three-element  Elastic-Viscoelastic-Plastic  model  is  superior 
to  those  made  using  both  the  standard  formulas  and  two-element  Elastic-Viscoelastic  Model. 

VII.  CONCLUSIONS. 

The  mechanical  behavior  of  materials  subjected  to  impact  loads  can  be  predicted  by 
using  relatively  simple  mechanical  models  to  describe  their  behavior.  The  model  parameters  can 
be  obtained  from  a true  stress-strain  diagram.  In  addition,  the  methodology  used  in  the  body  of 
this  report  can  be  used  to  predict  accurately  the  structural  integrity  of  munition  components. 
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GLOSSARY 


a stress 

e true  strain  = Cn(l  + eQ) 

eQ  engineering  strain 

E modulus  of  elasticity 

17  viscous  parameter 

6 plastic  modulus 

p density 

A shear  area 

I moment  of  inertia 

b thickness 

r radius 

v Poisson’s  ratio 

0 dynamic  shear  modulus  = 

2(1 +r>) 

dynamic  tensile  modulus 
w displacement 

t time 

A half  wave  length  = vQt 

vQ  impact  velocity 

3p  peak  acceleration 

T half  sine  wave  period 

F force 

M moment 
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radial  component  of  flat  plate 


circumferential  component  of  flat  plate 
initial  stress  wave  propagation  velocity  = 


radius  of  gyration  = 


wave  number 
pressure 

angular  frequency 

mass  X peak  acceleration 

mass 

factor  equal  to  (1  + »>/27r) 
amplitude  of  transverse  pulse 


group  velocity 


APPENDIX  A 


SOLUTION  OF  THE  DYNAMIC  EQUATION  FOR  THE  ELASTIC- 

VISCOELASTIC  MODEL  WHEN  EXCITED  BY  A HALF-CYCLE 

SINE  PULSE 


The  equation  of  motion  for  the  Elastic-Viscoelastic  Model  is: 


3e 


rearranging 


Vfo-  + (2E)e  = a 


be  + (2E)  € = o_ 
3t  n * V 


(1) 


or 


0e  + E'  _a_ 

at  rj  r] 


taking  the  Laplace  transform  of  this  equation 


sF(s)  + -jp  F(s)  - f(o)=“  Fc(s) 


(2) 


the  stress  generated  by  the  half  cycle  sine  pulse  is  m a (t)  = m ap  sin—  t or 

F(t)  - C0  sin(f)t 


where 


Cq  = m ap 


taking  the  Laplace  transform  of  the  force  equation  yields 


(3) 


IT 


Fo<*>  = Co“ 


s2  +(f  )2 


(4) 


substituting  equation  (4)  into  equation  (2)  yields: 

sF(s>  +irF(s)-f(o)  = jcofi 


s2+(f)2 
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rearranging  and  solving  for  F(s)  gives 


F(s)  = 


f + 0 

(0)  s2+  £)2 


E# 

s +| 


P,  ^ f(o)  + _S>  M 
F(s)  n'  TJ  UJ 


, E 

S+T 


E' 


(S+f-)(S2+f  )2) 


(5) 


The  time  dependent  strain  is  found  by  taking  the  inverse  transform  of  the  above  equation: 


e(t)  = 


jc-i 


1 

(s+-f  )(S2  +f)2 


If  we  neglect  the  residual  strain  on  the  system  prior  to  the  application  of  the  load  (eQ  = O)  the 
above  expression  simplifies  to 


e(t)  = ^f-f  X-1 


1 

(S+f-  )(s2+(f)2) 


expanding  the  fraction  in  bracket  by  partial 

eft)  (f)  J C“» 


performing  the  inverse  transform  yields: 


fractions  yields: 
Kj  K2 


k3 


s+-|j-  s+i(£-)  s-iff-) 
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where 


i = VT 


1 


substituting  these  constants  into  equation  (6)  yields 


e(t)  = 


+ 


rearrangement  along  with  using  the  relations;  cost  t = 1/2  + e-i(4)t) 

and 


yield  the  final  expression  for  the  strain: 


e(t)  = 


E'  , 
rf  sin1 

(?) 

It-* 

COS 

(f) 

I2  +1 

») 

2 

The  "dynamic”  modulus  is 


(7) 
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rearrangement  yields 


(t)2  + (If  »*(?)' 


y sin(r),-7  ““(?)*  +r  c"®' 


In  this  case,  the  burst  strength  was  to  be  measured  and  the  maximum  value  of  the  “dynamic  modulus’ 
was  of  interest.  This  occurs  when  sin  t = 1 or  t-j  , substituting  in  this  value  into  0 gives 


0 = 


2 . /'E’\2  1 


(?)  O 


f + ?e-(€>§ 


LT 


(8) 


when 


JL  E'  -/E'vr 

t > 1 7 and  e 


-d)f 


- l 


equation  (8)  can  be  simplified  to 

<t>  =(jr )ri 
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APPENDIX  B 


PROPAGATION  OF  STRESS  WAVE  IN  A FLAT  PLATE 


Figure  B shows  the  forces  acting  on  a small  element  of  the  plate  of  radial  dimension  AX. 
In  this  system  the  bending  moment  M is  balanced  by  shearing  forces  f acting  perpendicular  to  the 
plate  thickness. 


Figure  B 

From  Newton’s  second  law,  the  equation  of  motion  for  the  plate  is: 


. . 32w  df 

AxPA-p-  =f+Fx  Ax-f 


or 

a d2w  _ 3f 
p 3t2  3 x 

the  relationship  between  shearing  force  and  bending  moment  is 

dM 


(1) 


f = 


3x 
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where 


M = 01 


or 


f = 01 


3 

3X 


(2) 


where  w is  the  displacement,  x is  the  radial  component,  y is  the  circumferential  component  and  the 
strain  is  defined  as  equal  to 


3w  3w 

x=TZ-ey°!7 


The  relationship  between  the  radial  and  circumferential  components  are: 


v 


and  for  the  circular  flat  plate 


y = 2irx 


therefore 


3ey  _ 9ex 
3y  V by 


3ex  bX  _ v bex 
3 x 3y  27r  3 x 


substituting  the  preceding  equation  into  equation  (2)  yields 


±(^x  , - 3ex\ 
3x  \ 3x  2w  3x  / 


01 


+ 


JiAiffs. 

27T  / 3x  2 


or 


(3) 
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substituting  into  the  equation  of  motion  (1)  yields 


let 


w'W  Kt)S 


substituting  these  parameters  into  the  equation  of  motion  yield 


c02k2«2 


34w 

3x4 


The  solution  to  the  above  equation  is  w = Dcos(pt-ax)or  w = Dsin(pt-ax)  which  are  the  equations 
for  a transverse  pulse  where  a = p = — — and  A is  the  wave  length  of  the  pulse  and  c is  the  phase 
velocity  of  the  stress  wave. 

Substituting  either  one  of  the  above  equations  for  w into  the  equation  of  motion  gives, 


p2  = c02K2a4n2 


c = c 


o 


From  the  equation  of  continuity,  that  is,  the  displacement  of  the  disk  must  equal  the  total  elongation 
of  the  disk  or  vt  = ex.  Rearranging  this  equation  give  v=  ye  where  Von  Karman  showed  that 

X 

j is  equal  to  the  velocity  of  the  wave  propagation  in  the  plastic  region.  In  this  case  the  velocity 
of  wave  propagation  (Cg)  refers  to  the  velocity  with  which  a packet  of  waves  is  propagated,  the 
wavelengths  of  which  are  close  to  A.  This  is  related  to  the  phase  velocity  (c)  by  Cg=  2c. 

Finally  at  the  point  of  rupture 

0 " C°',,>  (irKO/ 
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APPENDIX  C 


CALCULATION  OF  THE  AISI  1018  LOW-CARBON  STEEL 
CIRCULAR  FLAT  PLATE 

From  the  stress-strain  diagram  (figure  6) 

A o = 65,000  - 4 1 ,700  = 23,300  psi 

eo’v  - ^eo’total  — eo,E^ 
eo,  total  = 0.30 
eo,E  = 002 

e0,v  = 0.30  - 0.02  = 0.28 


dgo 

at 


= cross  head  rate  X distance  between  jaws 


0-1  “/in/min 


in 


X 2 in  X 


t min 
60  sec 


= 0.00333  sec*1 


From  Calibration  Chart  (figure  5,  pad  31-200-085-350)  half  sine  wave  period,  r = 0.007  sec. 
From  test  data 


Acceleration  at  shear  = 145  g 
Disk  radius  = 1 in 
Disk  thickness  = 0.0153  in 

Handbook  Data 

E = 30,000,000  psi 
p = 0.284  lb/in^ 
v = 0.29 
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(1)  Calculation  of  dynamic  modulus: 

The  viscous  parameter  q,  was  calculated  by  substituting  the  values  given  in  figure  6 into 
the  equation  given  in  figure  2 of  the  body  of  this  report,  that  is: 


1? 


23,300 

.00333 


(1  + 0.28)  = 8,956,200 


lb-sec 


in 


2 


E _ 30,000,000  ...  .1 

1?  " 8,956,200  " 3,35  560 


» 3.1416  , 

r = moT  = 448-8secl 


E r 
t?  2 


= 3.35  = 0.0117 


since  >J1  an(j  e-0.0117  a ^ the  simplified  form  was  used, 


^=(7)  rj  = 448.8  X 8,956,200  = 4,018,500,000  psi 


However,  the  plate  was  in  shear  while  the  data  were  from  tensile  tests.  This  was  accounted  for  by 
using  the  following  relation: 


_ 4>E  4,019,500,000 

0S  " 2(1 +t>)  2(1+0.29) 


1,558,000,000  psi 


or 


224,352,000,000  psf  (lb/ft2) 

(2)  Calculation  of  stress  wave  parameters  where: 

- the  impact  velocity*  v = ap  ~ , in  this  caseap=  145  gs  acceleration 


or 


* See  Appendix  D for  response  characteristics  of  the  pulse 
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V 


145  g X 32.17  rt/soc2  X 


(0.007  sec) 
3.1416 


= 10.39  ft/sec 


- the  halfwave  length  A = vr 


A = 


10.39  ft/sec  X 0.007  sec  = 
- the  radius  of  gyration  K = 


.0727  ft 


1/2 


The  radius  of  gyration  is 


- the  propagation  velocity  c0 


= /224, 352, 000, 000 
\ 0.284  X 123 


1/2 


= 21,381  ft/sec 


(3)  Calculation  of  plate  stress  and  pressure  — from  Appendix  B. 


a = cnnv  I — — 
op  UrrKn 


= 21,381  X 490.75  X 10.39 


0.0727 


4 X 3.1416  X 0.188  X 1.023 
= 3,279,396  psf  or  22,774  psi  stress 
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The  pressure  P is  equal  to  the  total  load  (stress  X circumference)  divided  by  the  total  area  of  the 


plate,  P = 


o X 27rrb 
5 =° 

7rr 


for  a 0.0153-inch-thick  plate  having  a 1-inch  radius: 


697  psi  pressure. 


# 
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APPENDIX  D 


TUI  RESPONSE  CHARACTERISTICS  OF  THE  HALF-CYCLE  SINE  PULSE 

ACCELERATION 


a(t)  = ap 


t 


taking  the  Laplace  transform 


taking  the  inverse  transform 


v(t)  = vQt  + a 


= v0t  + a 


let 


v(t)  = “p(t)  ('-cos  ~)t) 

the  velocity  is  maximum  at  t = -~ 
having  vm  = ap(-J-). 
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US  Army  Science  & Technology  Center-Far  East  Office  1 
APO  San  Francisco  96328 

HQDA  DASG-RDZ  (SGRD-PL)  1 

WASH  DC  20314 

AMREL/MEB  i 

Wright-Patterson  AFB,  OH  45433 


Commander 

USEUCOM 

Attn:  ECJ%-0/LTC  James  H.  Alley 
APO,  NY  09128 
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